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Equal time, point to point correlation functions for spatially separated meson currents 
are calculated with respect to a variational construct for the ground state of QCD. Given 
such an ansatz we make no further approximations in the evaluation of the correlators. 
Our calculations for the vector, axial vector and scalar channels show qualitative agree- 
ment with the phenomenological predictions, whereas the pseudoscalar channel does 
not. However, the pseudoscalar correlator, when approximated by saturating with in- 
termediate one pion states agrees with results obtained from spectral density functions 
parameterised by pion decay constant and < —ipip > value obtained from chiral pertur- 
bation theory. We discuss this departure in the pseudoscalar channel , in context of the 
quark propagation in the vacuum. 

1. Introduction 

Quantum Chromodynamics (QCD) in the low energy sector is nonperturbative and 
the vacuum structure here is nontrivial.EJ The vacuum structure of QCD has been 
studied since quite^ome time both with quark condensates associated with chiral 
symmetry breakingp as well as with gluon condensates.™.! An interesting quantity to 
study with such a nontrivial structure of vacuum is the behaviour of current-current 



:erent spatial distances. This 
1 and studied through lattice 



correlators illustrating different physics involved at di ; 
has recently been emphasized in a review by Shuryak: 
simulations The basic point is that the correlators can be used to study the in- 
tcrquark interaction — its dependence on distance. In fact they complement bound 
state hadron properties in the same way that scattering phase shifts provide in- 
formation about the nucleon-nucleon force compicmcntary to that provided by the 
properties of the deuteron.13 

We have recently considered the structure of QCD vacuum with both quark and 
gluon condensates using a variational ansatz. We shall use here such an explicit 
construct of QCD vacuum obtained through energy minimisatiorfl to evaluate the 
meson correlators. 

We organise the paper as follows. In section 2 we recapitulate the results of 
Ref. 7. We then discuss in section 3 the quark propagation in our model of the 
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QCD vacuum. In section 4 we define and calculate meson correlation functions. In 
section 5 we quote the results. Section 6 is devoted to the study of the exceptional 
case of pseudoscalar correlator. Finally we discuss the results in section 7. 

We wish to make it clear that in this work we have calculated correlations of 
two currents each having the quantum numbers of the appropriate meson. The 
correlations are not between physical meson states.B 

2. QCD Vacuum with Quark and Gluon Condensates 

We have considered the vacuum structure in QCD using a variational approach 
with both quark and gluon condensates. Here we shall very briefly recapitulate the 
results of the same for the sake of completeness. The trial variational ansatz for the 
QCD vacuum is taken as 



\vac>=U G U F \0>, (1) 

obtained through the unitary operators U F and Ug for quarks and gluons respec- 
tively on the perturbative vacuum |0 >. 

For the quark sector, the unitary operator Up is of the form 

Up = exp(B F * - Bp), (2) 

with the quark antiquark pair creation operator Bp given by 

B F * = 



h(%y i(k)\a ■ jfc)cj(-fc) 



dk. (3) 



The operators S and c create a quark and antiquark respectively when operating on 
the perturbative vacuum. They satisfy the following quantum algebra in Coulomb 
gauge. 

[4 r (k),4 s (k') j } + = 6 rs 5^S(k- k') = (4) 

Further h(k) is a trial function associated with quark antiquark condensates. 

Clearly Eqs. (||) and (||) correspond to operator equations which create an ar- 
bitrary number of quark antiquark pairs. In fact Eq. (^) may be interpreted as an 
operator to create a Bose BCS state. We shall further assume h(k) to be spherically 
symmetric so that the non perturbative vacuum state will have the same symmetries 
as the perturbative vacuum i.e. zero momentum and angular momentum. 

Similar considerations are applied for constructing the gluon condensate func- 
tion. Thus for the gluon sector we have, 

Ug = exp (Bg^ — Bg) ; b g^\ J f(k)a a t (k) f a a t (-k^ dk. (5) 

Here f(k) is a trial function associated with gluon condensates and a a i{k) the 
transverse gluon held creation operators. 
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Clearly such a structure for the vacuum eventually reduces to a Bogoliubov trans- 
formation for the operators. One can then calculate the energy density functional 
given as 

e = F(h(k)J(k)). (6) 

The condensate functions f(k) and h(k) are to be determined such that the energy 
density eo is a minimum. Since the functions cannot be determined analytically 
through functional minimisation except for a few simple cases,! we choose the al- 
ternative approach of parametcrising the condensate functions as ( with k — \k\), 

A' 

tan2fr(fc)= (ei?2fc2 _ 1)1/2 - (7) 

This corresponds to taking a Gaussian distribution for the perturbative quarks in 
the nonperturbative vacuum. 

Similarly, for the function f(k) describing the gluon condensates we take the 
ansatz sinh/(fc) = Ae~ Bk I" 1 . 

Further one could relate the quark condensate function to the wave function of 
pion as a quark antiquark bound statei and hence to the decay constant of pion. 
It is clear from Eq. (ffl) and the ansatz for the gluon condensate function that even 
for k — 0, f(k) and h{k) are non zero. Hence there is a finite probability to find 
low momentum states in the vacuum. 

In Ref. 7 the energy density is minimised with respect to the condensate param- 
eters subjected to the constraints that the pion decay constant f^ and thfi gluon 
condensate value — < G° , ,,,G atlu > of Shifman Vainshtcin and ZhakarovEJ come 
out as the experimental value of 93 MeV and 0.012 GeV 4 respectively. 

The results of such a minimisation showed the instability of the perturbative 
vacuum to formation of quark antiquark as well as gluon condensates when the 
coupling became greater than 0.6. Further the charge radius for the pion comes out 
correctly (R c h — 0.65 fm) for a s = 1.28. The corresponding values of A and R of 
Eq. (0) are calculated to be A' min ~ 1 and R ~ 0.96 fm. 

It should be pointed out that local colour neutrality is lost in the vacuum struc- 
ture we have due to the generation of mass like terms which is a limitation of our 
approach. 

With the structure of QCD vacuum thus fixed from pionic properties and SVZ 
value we consider quark propagation in the next section. 

3. Quark Propagation in the Vacuum 

In the calculation of correlators, quark propagators enter in a direct manner 
and hence it is instructive to study aspects of the interacting propagator in some 
detaiLEil The reason for doing this is two folds. We wish to know how it differs from 
a free massive propagator i.e. how good it is to have a "constituent quark" picture 
and further to compare and contrast with other approaches such as instanton liquid 
model or vacuum dominance model based on operator product expansion. 
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The equal time interacting quark Feynman propagator in the condensate vacuum 
is given as 

s«p($) = (^[i>i(z),$m]Y (8) 



In our model this reduces to 



11 / ik. 



S(x) = -— -r e lk:x dk sin2h(k)~tf -It) cos2h(k) (9) 
2 J L J 

- ' ' " ' 1 1 e^^--^^I(x), (10 ) 



2tt 2 x a (2tt) 3 / 2 2i? 3 (2tt) 2 x 



where 



DC: 



/LH = / ( COSfcr " 1+(1 fc ! e -fl^)l/2 dfc ' ( U ) 



(I 



and x = fc = |fc|. 

Clearly, the free massless propagator is given by Sq(x) = —^-^r which can be 
derived independently or from ( Eq. ([h])) in the limit R : — » oo i.e. in the limit that 
the condensate functions vanish. 

It may be useful to note that the free massive propagator is given as 

1 771 2 

S (m q ,x*) = 9 [Ki(m q x) - i 7- xK 2 {m q x)] , (12) 

\2tt) 1 x 

where K\{m q x) and Ki(m q x) are the first and second order Bessel functions re- 
spectively. 

In Fig. ^we plot the two components Tr S(x) and Tr (^■x)S(x) of the propagator 
for massless interacting quarks given by Eq. (|lo| ) corresponding to the chirality 
flip and non-flip components considered by Shuryak and VerbaarschotJid The first 
trace is normalised to the short distance limit of the massive free quark propagator 
Tr So(m q ,x)/m q which from (Eq. (|l2|)) is l/(7r 2 a: 2 ). The second trace is normalised 
to the free quark propagator which is Tr (7 • x)Sq(x) — 2 i /(ir 2 x 3 ). 

To compare with the constituent quark models with an effective constituent 
mass, we have also plotted the behaviour of free massive quark propagator with 
masses 100 MeV, 200 MeV and 300MeV. In the chirality flip part, the propagator 
in the condensate medium starts from zero, consistent with zero quark mass at 
small distances, attains a maximum value of about 200 MeV at a distance of about 
1.3 fm and then falls off gradually. Further the interacting propagator overshoots 
the massive propagators after about 0.8 fm. These features are qualitatively similar 
to that of the instanton liquid model of QCD vacuum considered in Ref. 11. 

In the chirality non flip part, the interacting propagator starts from 1, again 
consistent with the behaviour expected from asymptotic freedom, but at larger sep- 
aration it falls rather slowly indicative of an effective mass of the order of 150 MeV. 
Again these features are similar to that of Ref. 11, though quantitatively there are 
differences. 
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7T 2 x 2 Tr S(x) [in MeV] tt 2 x 3 /2 Tr[y . x S(x)] 




0.5 1 1.5 0.5 1 1.5 2 



x (in fm) x (in fm) 

Fig. 1. The two components of the quark propagator. Tr(S) (a) and Tr[(7 • x)S] (b) versus the 
distance x(in fm). The normalisation, indicated in the figure, have been explained in the text. 
The solid line corresponds to massless quark interacting propagator S(x). The three lines, short 
dashed, dot-short dashed and long dashed correspond to a massive free propagator with a mass of 
100, 200 and 300 MeV, respectively. 

It may be amusing to consider the leading behaviour of the propagator as x — •> 0. 
In this limit, the first term of Eq. (|^) is given by 



Ti = J dksin2h(k) + 0(x 2 ) (13) 

= ^<~^>+0(x 2 ), (14) 

where < — ipij> >=< uu + dd > for two flavours. Similarly the leading contribution 
from the second term of Eq. (^|) is 

T^-^ + 0(x). (15) 
Thus in the small x limit the interacting propagator of Eq. (^) reduces to 

m=-^^+±<-w> as) 



Quark Propagator and Meson Correlators in the QCD Vacuum 6 

which is exactly the result of the vacuum dominance model0 based on operator 
product expansion. 

4. Meson Correlation Functions 

Consider a generic meson current of the form 

J(x)=rPi(x)T a0 ^ (x) (17) 

where a; is a four vector; a and /3 are spinor indices; i and j are flavour indices; T is 
a 4 x 4 matrix (1,75, 7m or 7m 7s) 

Because of the homogeneity of the vacuum we define the conjugate current to 
the above at the origin as, 

j(o) = ^ A (o)r' A ^4(o) (is) 

with r = 7oT^7o 

In general, the meson correlation function is defined as, 

R(x) =< TJ(x)J(0) > vac . (19) 

From now on we assume that expectation values are always with respect to the 
nonperturbative vacuum of our model, hence we drop the subscript vac. 
Hence with Eqs. @, © and (|l|) we have 

R( X ) = r a0 r' xs < r$,(aO^(aO$(o)$(o) > . (20) 

This reduces to the identity 

R(x) = r a pf xs < ZVj(*)$(0) X Tiji(x)4(0) > . (21) 

The above definition of R(x) is exact since the four point function does not con- 
tribute. In fact, in the evaluation of Eq. ( |20| ) we shall have a sum of two terms. 
The first is equivalent to the product of two point functions which is Eq. (|l]) . The 
second term arises from contraction of operators at the same spatial point, related 
to disconnected diagrams and thus can be discarded. 

In Eq. (|^j) the first term can be identified as the interacting quark propagator 

S(x) =< T^(x)$ j (0) > 

It can be shown using the CPT invariance of the vacuumEl that the second term 
is given as 



<T^{xW(Q)> = -755(0075 



(22) 
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Hence the correlation function of Eq. ( pp| ) becomes 

R(x) = -Tr \s(x)T S(-x)T~\ . (23) 
Similarly the correlator for massless noninteracting quarks can be given as 

Ro(x) = -Tr [s (s)r'So(-aOr] . (24) 

Our task is now to evaluate the expression ( p3[ ) with the ansatz for QCD vacuum 
as given in Eq. (^). Since we have evaluated the interacting quark propagator at 
equal time we also consider the correlation functions (Eqs. (p3|)-(p^)) at equal time. 
This implies that the four vector x in the above equations reduces to the three 
vector x 

Having obtained the propagators in the earlier section, we can calculate the 
correlation function, Eq. ( p3| ) for a generic current of the form as in Eqs. ([l7]) and 

For convenience, we will consider the ratio of the physical correlation function 
to that of massless noninteracting quarks and obtain (with x = \x\) 



R(x) 
Ro(x) 



1 + -x 2 I(x) 



2 6 2i 2 x 2 Tr 



r r 



8R e 



% x iTr [7 i rvr] ' 



(25) 



which is then evaluated in different channels with the corresponding Dirac structure 
for the currents. 

Table 1. Meson currents and correlation functions 



CHANNEL 


CURRENT 


PARTICLE 
(J P ,MASS in MeV) 




COI 


IRELAT 
' R(x) 
Ro(x) 


OR 


a 


Pseudoscalar 


JP = u'jsd 


7r°(Cr,135) 


[n 


-\x 2 I(x)] 2 + l 


x 6 

W e 


-x 2 /R 2 


Scalar 


J s = ud 


none(0 + ) 


[n 




x 6 

w e 


-x 2 /R 2 


Vector 


Jfi = u^^d 




[n 




-x 2 /R 2 


Axial 


Jl = ""7^75 d 


Ai(l+,1100) 


[n 


-i x *I( x )] 2 -^e 


-x 2 /R 2 
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5. Results 

We have studied the above ratio of correlators for four channels. In each channel 
we associate the current with a physical meson having quantum numbers identical 
to that of the current. The results are shown in Table [j] and in Fig |[ 




0.5 1 1.5 2 0.5 1 1.5 2 

x (in fm) x (in fm) 



Fig. 2. The ratio of the meson correlation functions in QCD vacuum to the correlation functions 

R(x) , „ 

for noninteracting masslcss quarks, , Plotted vs. distance x (in fm) 

Ro(x) 

We may notice some general features and relationships among the correlators. 
The pseudoscalar correlator is always greater than the scalar correlator and vector 
correlator is greater than the axial vector correlator. We may emphasize here that 
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these relations are rather general in the sense that they do not depend on the explicit 
form of the condensate function and arise due to the different Dirac structure of 
the currents which is reflected in the generic expression for the correlation functions 
as in Eq. ( |25| ) . The behaviour of each channel is consistent with that predicted by 
phenomenology except in the pseudoscalar case where the ratio does not go as high 
as expected from phenomenology. We examine this in the next section. 

6. Pseudoscalar Channel 

The explicit evaluation of the pseudoscalar correlator gives, using Eq. (p5|) 



R(x) 
Rq(x) 



1 + \x 2 I{x) 



which may also be read off from column 4 of Table [l| This is plotted as a function of 
x in (Fig. |^). As may be seenfrom (Fig. ||) this ratio has a maximum of ~ 1.2 at x ~ 
1.3 fm. PhenomenologicallyO the peak is at ~ 100 at x ~ 0.5. In order to compare 
our results with other calculations we evaluate the same correlator approximately 
by saturating intermediate states with one pion states. 

With our definition of the correlation function (Eq. (|l9|)) we have 

R{x) = < TJ p (x)J p (0) > (27) 
- ~ (< J p {x)J p (0) > + < J p (0)J p (x) >) (28) 

which comes from the definition of time ordered product for equal time algebra. 

We now insert a complete set of intermediate states between the two currents 
but retain only the one pion state in the sum for the four point function. Thus, 



R(x) = ~ J {< vac | J p (x) | n a (p) >< n a (p) | J p (0) | vac > 

+ < vac | J p (0) | n a (p) >< n a (p) | J p (x) | vac >)dp (29) 

Using translational invariance and the fact that for the pseudoscalar current 
J p = U75C? and J p — —J p , the correlator may be written as 



R(x) = ^J<vac\ J p (0) I TT a (p) >< TT a (p) I J p (0) I vac > (e lp - s + e~ lp - s ) dp 

(30) 

We may evaluate the above matrix element using the definition of the pion decay 



constant given asEJ 



where J pa — [ip^ tl 'y 5 T a ip] is the axial current. 
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It can be shownEl that the divergence of the axial current gives the pseudoscalar 
current 

d„ [V^V^] =2im, [^ 7 5 r a ip] (32) 

where m q is the current quark mass. Thus taking divergence of both sides of Eq. 
© and using Eq. (§|) we get, 

— f 711 ^ 

2 m q < vac | zJ-(x) | n"(p) >= ^j^,^ (33) 

where we have used pL = m\ . 

In an earlier papera within our vacuum model and using the fact that pion is an 
approximate Goldstone mode it was demonstrated that saturating with pion states, 
gives the familiar current algebra result 

ml = < -fa, > (34) 

Jit 



With this result we eliminate quark mass m q in Eq. (|33|) in favour of the quark 
condensate to get the relation 

<"°'^<°)'^>° 2(2 T )»4p,)^ < T > <35> 
Using this, the expression for the pseudoscalar correlator now becomes 

- eh Mr^) 2 / wrk^ {e "' + ,HM > * (36) 



The above integral can be evaluated using the standard integraS 
r p (p 2 + 2 Y- 1 l 2 sm{ap)dp = 4= ( — 

Jo v 71 " V a 



' 1 cos(vk)T(v + i)Jf„ + i(a/3) 



for a > 0, Re/3 > and in the limit v — » 0. We then finally get for the correlator 
(using saturation of pion states and with x = \x\) 



= is? < 37 > 

The correlator for free massless quarks as calculated in the earlier section for pseu- 
doscalar is 

R (x) = -±-g (38) 



Hence the ratio is 

R(x) 7T 2 / < —^ip > 



Ro(x) 16 V U 



x m^Kxim^x) (39) 
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Fig. 3 (a) 




0.5 1 1.5 2 

x (in fm) 



Fig. 3. The pseudo-scalar correlator plotted for < 
that of Shuryak < -ipip >= (307.4 MeV) 3 in (b) 



Fig. 3 (b) 




0.5 1 1.5 2 

x (in fm) 



value of < —ipip >= (190 MeV) 3 in (a) and 



We have^pkitted in Fig. || this ratio for our value of < —ipip > and that used 
by ShuryakBlia Note that our value of < —Tpiji > is an output of the variational 
calculation consistent with low energy hadronic properties^ We thus observe that 
the approximate calculation of the pseudoscalar correlator due to saturation with 
one pion states (Fig. ||(a)) yields higher values (~ 15 times more) as compared 
to the calculations without saturation as aiL approximation (Fig. |J) . Thus the 
fermionic condensate model for QCD vacuumQ does not give as high values for the 
pseudoscalar correlator as required by phenomcnological results. The value we have 
used for < -^V >- (190 MeV) 3 is smaller than Shuryak's value of (307.4 MeV) 3 b 
which appears in the parameterisation of the physical spectral density through the 
coupling constants With his value of < — -ipip > the ratio R(x)/Rq(x) is shown in 
(Fig. ||(b)) which agrees with phenomenology. 

7. Summary and Discussions 

We have evaluated the mesonic correlators in this paper using a variational con- 
struct for the QCD vacuum. Except for the pseudoscalar channel the results show 
qualitative agreement with phenomenological results a We have also shown that the 
quark propagation in our construct for the QCD vacuum is almost identical to that 
in the instanton modeL 

Following ShuryakJjO we also see that using current algebra approach the pseu- 
doscalar correlator rises sharply with spatial separation. Let us recall that the cur- 

6 Our definition of the condensate value differs from the standard one by a factor of 2 1 / 2 
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rent algebra result also follows from the approximation of saturating by one pion 
states in the normalisation of the pion state.a 

It might appear that by suitably changing the value of < —ipip > in our calcula- 
tions without saturation one might be able to reproduce all the phenomenological 
results. Actually we find that it is not so. In fact, it adversely affects the correla- 
tors in the other channel which can be seen in the expressions given in column 4 of 
Table |l|. Also the results of our calculations are not very sensitive to the parameter 
R. 

In view of these findings, it is not clear whether saturation of intermediate 
states by one pion states only in the evaluation of the correlator, is sufficiently well 
justified. We therefore think that a unified treatment of correlation functions in all 
the channels is still not available. 
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